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Abstract: Abraham's force seems to have been observed in very low frequency experiments, 
but the existence of an Abraham's force exerted by light is still to be proved. In fact, the only 
experiments performed with light have measured the radiation pressure exerted at the interface 
between different dielectric media, and their value for inferring conclusions about the existence 
of Abraham's force is at best opinable. Twenty years ago Brevik proposed an experiment for 
checking what forces light exert in the interior of a transparent medium. It is argued that his 
proposal is feasible today with a detection method now at hand, thanks to the great improvements 
achieved in manufacturing optical fibres with the required physical properties. 

1. - Introduction 

Let us figure out a physicist who has engaged in even a cursory survey of the huge 
amount of theoretical literature accumulated in more than one century about the forces 
that electromagnetism supposedly produces in dielectric media: he has patiently compared 
the divergent outcomes that eminent theoreticians have arrived at; he has doubtfully pon- 
dered the relative advantages of the microscopic versus the macroscopic approach; he has 
followed the alternating fortunes of electrostriction and magnetostriction 1 ); he has won- 
dered in disbelief why the supposedly universal tool of quantum mechanics has been so 
inept, in a time span of seventy years, in shedding some light on this issue. While chal- 
lenged by such a confusing state of affairs, he has even been summoned to surrender, and 
recognize the evident vanity of the whole effort 2 ' . 

No wonder whether this physicist, like the simple Simon of the stinging tale [11] 
written by J.L. Synge and entitled "On the present status of the electromagnetic energy 
tensor", might feel that his very perception of an objective reality is at stake, and that, as 
an antidote, he had better giving a look to whatever relevant information experimentalists 
may have gathered on these forces, in the hope that the latter have turned out to behave 
like classical, macroscopic forces usually do, and that they have not been found to depend 
on the mental act of splitting the overall energy tensor in two. 



Dedicated to Prof. Hans-Jiirgen Treder on his 70th birthday. 

+. Contributed to: .From Newton to Einstein (A Festschrift in Honour of the 70th Birthday of 
Hans-Jiirgen Treder), W. Schroder Editor, Bremen: Science Edition 1998. 

^ Accounted for already in the seminal works by Helmholtz [1],[2], they have no place in the 
relativistic proposals set forth by Minkowski [3] and Abraham [4], [5], but are present (with a 
different shape) e.g. in the paper by Einstein and Laub [6]. 

2 ) For instance, in Refs. [7]- [10] it is argued, under natural assumptions like those introduced 
in Ref. [9], that the way one splits the overall macroscopic energy tensor into a matter term and 
a radiation term is just a convention, merely dictated by reasons of convenience, and maybe of 
elegance. 
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2. The experimental evidence for Abraham's force 



Given the huge amount of theoretical literature, our natural philosopher might expect 
to be confronted with an equally fair production on the experimental side, but to his 
dismay, he will discover that all the papers dealing with experimental findings can be 
easily contained in a small drawer 3 ). 

Clear evidence for the Abraham force 4 ) was retrieved with two kinds of experiments, 
both availing of low frequency electromagnetic fields [15]-[17]. One experiment performed 
by Walker and Lahoz will be summarized here, since it seems [12] to have provided the 
most unambiguous outcome up to date. 

To do so, let us first specialize the expression of Abraham's generalized force density 
given in the Appendix for the case when a global coordinate system exist in which the 
components of the metric are: 

gik=rjik = diag(l,l,l,-l), (1) 
and the components of the four-velocity of matter are 

u 1 = u 2 = u 3 = 0, u 4 = 1. (2) 
Let us assume also the vanishing of the four-current s l ; then eq. (A33) reduces to: 

f P = -^[F 4A F/] 4 , f 4 = 0, (3) 

where e and fi stand respectively for the dielectric constant and for the permeability of a 
homogeneous and isotropic medium at rest. Greek indices run from 1 to 3 and label the 
spatial components, while a comma signals ordinary differentiation. Then in m.k.s. units 
Abraham's ordinary force density takes the well known expression: 



f = 



n 2 - 1 
c 2 di 



(EAH), (4) 



where n is the refractive index of the medium, and c is the velocity of light in vacuo. 

In the experiment [16] of Walker and Lahoz a disk of barium titanate, whose relative 
dielectric constant was near to 4000, had a small cylindrical hole pierced at its center, 
and both the inner and the outer cylindrical surfaces were coated with metallic films. 
By sending charges on these electrodes a radial electric field could be set up within the 
disk. The latter was suspended to a tungsten wire, and the torsional pendulum obtained 



3 ) See, e.g., the accurate report [12] published by Brevik in 1979; to our knowledge, just one 
further finding related to our issue, about the angular momentum of radiation in dielectrics [13], 
has occurred in the meantime. 

4 ) A simple derivation of Abraham's energy tensor, found by W. Gordon [14] in 1923, is outlined 
in the Appendix. 
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in this way was subjected to a uniform magnetic field of constant value, parallel to the 
suspending wire. When the radial electric field was allowed to vary according to a sine law 
at the resonance frequency of the pendulum, steady oscillations were observed to set in, and 
their amplitude was measured with an optical lever. From the known parameters of the 
experimental device a value for f exp . was inferred, in good agreement with the theoretical 
prediction of eq. [4], as far as can be achieved in a delicate experiment of this sort. 

3. Reviving a proposal by Brevik 

It would be interesting to learn whether light exerts Abraham's force on a transparent 
medium crossed by it. Since in eq. (4) the product E A H is differentiated with respect 
to time, a large value of the instantaneous force density is predicted in this case, but the 
direct observation of a mechanical force at the frequency of light is beyond reach also for 
present day experimenters. Brevik [18] made some proposals for detecting the Abraham's 
force of light in an indirect way, and one of them will be reconsidered here. An order 
of magnitude estimate of several quantities intervening in the experiment is needed for 
assessing whether the latter is actually feasible. 

Imagine a single-mode optical fibre 5 \ whose length is L = 1km. The mathematical 
description of the eigenmodes of the electromagnetic field propagating in such a fibre is 
rather involved [21], but if we content ourselves with an order of magnitude evaluation, 
the actual field for the so called transverse electric mode can be replaced with the field of 
a plane wave propagating along the axis of the fibre, whose intensity is constant within 
the fibre core, and vanishes abruptly at the core radius. Let the axis of the fibre coincide 
with the x axis of a Cartesian co-ordinate system, and the electromagnetic plane wave be 
a linearly polarized one 6 ) with the nonvanishing components: 

E y = Eocos(kix — uJit)cos(k2X — a^t), (5) 
B z = ^eoenonE cos(kix — uii)cos(k2X — u^t), (6) 

where ui\ <C and 

k\ ki n 



5 ) For light with a wavelength in vacuo A = 1.3/Um, as will be chosen here, a single-mode fibre 
is typically constituted by a cylindrical doped glass core with a radius T\ = 5fJ,m, surrounded by 
a cylindrical shell of pure silica glass with a slightly lower refractive index of radius T2 = 60/Um, 
and by a further cylindrical coating made of acrylic resin of still lower refractivity, up to an outer 
radius = 125/Um. For an account of the way optical fibres behave see e.g. Refs. [19] and [20]. 

6 ) Polarization preserving optical fibres are a recent achievement; in usual fibres the polarization 
is elliptical with irregular variations of the parameters of the ellipse as one proceeds along the 
axis. This circumstance will not affect our order of magnitude estimate. 
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i.e. we imagine sending 7 ) into the fibre linearly polarized light whose amplitude is modu- 
lated with the low angular frequency Let us set \x = 1, as appropriate for glass. Then 
Abraham's force density, defined by eq. (4), is directed along the axis of the fibre and 
takes the value: 

n 2 — 1 d 

f = n— e £ 2 x- {[1 + cos2{k x x - vit)] [1 + cos2(k 2 x - u 2 t)]} . (8) 

The right-hand side of eq. (8) contains a nonvanishing component at low frequency that 
is potentially observable; it reads: 

n 2 - 1 

flow = n — e E <jj 1 sin2(k 1 x - w\t). (9) 

2c 

The average intensity of the radiation described by eqs. (5) and (6) is I = |e ci? 2 . Let us 
assume for Eq the value 

E = 10 6 ^; 

m 

then the previously quoted size of the core radius allows the fibre to carry the power of 
~ 60mW, an acceptable choice. Imagine now that the fibre is wound, together with a 
second one, on a cylinder with radius R = 10cm; this move looks possible without altering 
in an appreciable way the fields prevailing in the core. The cylinder is suspended to such 
a tungsten wire as to constitute a torsion pendulum whose resonance angular frequency u> 
is 



u = 2ui = 2ixs 

Due to the low frequency component of Abraham's force the first fibre will exert on the 
torsion pendulum the torque 8 ^, say: 



M = irr\Li low R = M Q sinut, (10) 

where, due to the previous choices, and since in the core n = 1.44 at the wavelength chosen 
for the light: 



M = 0.63 x IQ-^Nm. 
The deflection angle d of the torsion pendulum shall fulfil the equation 

H + d-d + kd = M sinujt, (11) 

7 ) A fibre laser [20] can be directly joined to the fibre without interposition of reflecting surfaces. 
Eqs. (5) and (6) do not account for losses; this is a reasonable approximation, since the attenuation 
of the fibre considered above can be as low as 0.3 dB, the limiting value due to Rayleigh scattering. 

8 ) If both ends of the fibre are left loose from the torsion pendulum and rigidly clamped, surface 
forces occurring at the extremities will not contribute to the torque. 
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where /, d and k are respectively the momentum of inertia, the damping parameter and the 
elastic constant. For very small damping the resonance angular frequency uj is practically 
equal to the natural frequency of the undamped pendulum, i.e. uj = y/k/I. The standing 
oscillation at that frequency, according to eq. (11), obeys the law: 

= -Ml cosut . (12) 
ud 

We can estimate the damping parameter d by reminding that the time constant r of the 
damped oscillations is related to d by the equation r = 21 /d. Since a reasonable value for 
the time constant is r ss 10 3 s, while / « 2.5 x 10~ s kg x m 2 , the standing oscillation at 
resonance will occur with the tangential velocity: 

v = R$ = — '^sinujt, 
d 



where one can set 

d 



1.26 x lO^ms -1 . 



We propose availing of the second fibre to detect this velocity. We shall disregard 
the fact that the fibre, whose length can be again L = lkm, is wound on a cylinder, and 
reason as if it had been uncoiled. We shall assume that when this fibre is at rest its core, 
parallel to, say, the x axis, is travelled by the unmodulated, linearly polarized light 9 ) 
whose components, written more relativistico, are: 

F 2 4 = E' cos(k 2 x - uj 2 t), (13) 
F12 = y/ejIE' cos(k2X - u 2 t). (14) 

We imagine however that the velocity field has followed the fibre in the deformation due 
to the unwinding action; given the smallness of v, we can approximate the four-velocity of 
the uncoiled fibre as: 

0,0,1), (15) 

where 

M R . M R . 

p = — smut ~ — sinlzl k\x — uj\t) J, (lb) 

cd cd 

because k\x is a negligible term also for x = L. Since the fibre is subjected to this velocity 
field, eqs. (13) and (14) no longer provide an exact solution to the Maxwell's equations: 
the constitutive relation has changed, however slightly, in keeping with eq. (A14). Due to 
the smallness of the change, calculating a first order perturbation may suffice. If 5F rs is 



9 ) F or reasons apparent in the sequel, the second fibre needs to be polarization preserving. 
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the change undergone by the unperturbed F rs , the corresponding change in H rs will be 
given by: 

11 1 

5H rs = -5F rs + e -^^(5l5F 4s -5i5F 4r ) + e -^^[3(5 r 1 F 4s -5tF 4r +5lF sl -5 s 4 F rl ). (17) 

Hence in the present case the first order correction to one unperturbed set of Maxwell's 
equations reads: 

-8F x f p + e5F x j = e -^—± [((3F 4X ) A - (/3F 1A ), 4 )] (18) 

for r = A, and 

e5F% = 0, (19) 
for r = 4, while the correction to the other unperturbed set is: 

$F [ikM = 0. (20) 

Eqs. (18)-(20) display in 8F ik the same form as the unperturbed equations do in with 
the exception of the component with A = 2 of eq. (18), that reads: 

5F 2 \ + e^5F 24 = 2(e^ - 1)((3F 42 ) A (21) 

since in our case: 

(/3F 42 ) ;1 = -(/3F 12 ), 4 . 
Therefore only 5F 21 and 5F 24 need to be nonvanishing. Let us set: 

C = (2fci + k 2 )x - (2wi + u 2 )t, 
t] = (2ki - k 2 )x - {2uj\ - u> 2 )t, 

and 

then a physically appropriate particular solution is given by: 

5Fi2 = C {sin( + sint] + (2k± + k 2 )xcosQ + (2k± — k 2 )xcosrj} , (22) 
C 

5F 42 = — {(2ki + k 2 )xcos( + (2&i - k 2 )xcosr]} . (23) 



We are interested in the behaviour of SF^ near the end of the second fibre, therefore we 
can retain only the terms that grow linearly with x. Since k± <C /c2, we eventually write: 
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IvI RE 1 x 

SF12 ~ 27r(e / u — 1) —^- — sinut sin(k 2 x — u 2 t), (24) 

cd A2 

c-7-, c efi — 1 M RE' x . ... . 

0^24 ~ 27T — : — SlUUjt Sin\K2X — UJ2I) , (25) 

^/e/Z ca A 2 

where A2 is the wavelength of light in vitro. We remark that SFik, at variance with 
the unperturbed F^, exhibits a "sin" rather than a "cos" dependence on the argument 
k2X — U2t. The detection of the signal could be done by adding 10 ) at the end of the second 
fibre the field with the nonvanishing components: 



F12 = v ^ejIE' sin(k2X - u 2 t). (26) 

F 2 4 = E' sin(k 2 x - u 2 t). (27) 

Then the square of the electric field E y has a component whose low frequency average, 
that can be perceived by a photo detector, reads: 

{E 2 y ) av . = (E' ) 2 [l - 2. x 10- 10 sinu4 (28) 

Although quite small, when compared with the background, the signal caused by the 
component of (Ey) av _ at the resonance frequency of the torsion pendulum seems within 
reach of the present day techniques of digital phase-sensitive detection. Its observation 
would provide evidence that Abraham's force is indeed exerted by light inside a dielectric; 
in fact, electrostriction [7], [12] in glass should produce a signal that is smaller in magnitude, 
and opposite in sign with respect to the one calculated above. 
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10 - ) Fibre couplers just built with this scope in mind are commonly available. Of course the 
whole procedure of detection is only possible if the second fibre is polarization preserving, and if 
the coherence length of the light sent in the second fibre is longer than the fibre itself. The latter 
condition ultimately means [20] a request about the phase and the amplitude stability of the laser 
source that can be met with by present day lasers. 
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Appendix: Gordon's derivation of Abraham's tensor 



The general relativistic argument conceived by W. Gordon has been retrieved in the 
limbo of the forgotten papers [22]. It runs as follows. In a four-dimensional differentiable 
manifold let us consider a contravariant skew tensor density H lA: , a covariant skew tensor 
Fik, and write the naturally invariant equations: 

Hj = s\ (Al) 



F [iKm ] = 0. (A2) 

The vector density s l , defined by the left-hand side of eq. (Al), represents the electric 
four-current density, while the comma signals ordinary differentiation, and we have set 
F[ik,m] = l(Fik,m + Fkm,i + F mi , k ). Equations (Al) and (A2) express Maxwell's equations 
in general curvilinear co-ordinates; they need to be complemented with the constitutive 
relation of electromagnetism, a tensor equation that uniquely defines e.g. H lk in terms of 
F^ and of whatever additional fields may be needed for specifying the properties of the 
electromagnetic medium. When the assumed dependence of H lfc on F^ is algebraic and 
linear, the constitutive relation reads: 

H lk = iX lfc ^F mn , (A3) 

and the electromagnetic properties of the medium are summarized by the four-index tensor 
density X ifcmn . In the case of vacuum one writes: 

XfjT = vWV n - 9 tn 9 krn ); (A4) 

the constitutive relation entails only the metric tensor g ik and g = —det(gik) in the way 
known from general relativity. In this case Maxwell's equations are usually written as: 

F$ = s\ (A5) 

F [ik , m ] = 0, (A6) 

in terms of a skew tensor F ik and of the contravariant tensor density F lk = ^/gg irn g kn F mn . 
A general skew tensor Fik can always be written as the sum of the curl of a potential and 
of the dual to the curl of an "antipotential" : 

Fik = <Pk,i - <Pi,k + e ik mn {^ n ,m ~ 1pm,n), (A7) 

where e ik mn is the tensor obtained from the Ricci-Levi Civita symbol e lkmn in the usual 
way. Starting from the Lagrangian density 

L = \¥* k F lk - sVi (A8) 

for the general relativistic vacuum, both sets of Maxwell's equations can be derived through 
the Hamilton principle, by asking that the variations of A = J LdS (dS = dx 1 dx 2 dx 3 dx 4 ) 
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with respect to <pi and to ipi separately vanish [23] . The well known form of the energy ten- 
sor density Ti k for the electromagnetic field in vacuo is obtained through the method [24] 
inaugurated by Hilbert, i.e. by carrying out the Hamiltonian derivative of the Lagrangian 
density (A8) with respect to g lk : 



T lk = 2^ = F^F fcn - \g ik F mn F mn . (A9) 



SL 

'Jg~i 

The derivation of the energy tensor must be performed by keeping into account the homo- 
geneous field equations obtained above; they dictate that F ik is the curl of a four-vector 
<fi and does not contain the metric. 

The constitutive relation for a linear, nondispersive medium will have in general the 
form of eq. (A3). Gordon deals with a medium that in addition is homogeneous and 
isotropic, when considered at rest; its constitutive relation can be written in a simple 
form, due to Minkowski [3] , that can be extended without change to general relativity. Let 

i = 

be the four- velocity of matter, for which UiU 1 = — 1. One defines the four- vectors 

Fi = F lk u k , ^ = H ik u k , (All) 

where H ik = (1/ y/g)g ip g kq HP q is the covariant tensor associated with H lh . Then the above 
mentioned constitutive relation simply reads: 

Hi = eF,, (A12) 

UiF km + u k F mi + u m F ik = n{uiH km + u k H mi + u m H ik ), (A13) 

where e and fx are the dielectric constant and the permeability, and do not depend on the 
chosen event. These eight equations, that entail two identities, are equivalent [25] to the 
six equations: 

fxH ik = F lk + {en - l){u l F k - u k F l ). {A1A) 
Since the right-hand side of eq. (A14) can be rewritten as 

F rs {g ir 9 ks - (e/j - 1) (u i u r g ka + u k u s g* r ) } , 

and, due to the antisymmetry of F rs , one can freely add the term (eii~l) 2 u l u k u r u s within 
the curly brackets, the constitutive relation eventually comes to read: 

nW k = [g ir - (en - l)u l u r ] [g ks - (e/x - l)u k u s ] F rs . (Alb) 
The "effective metric tensor" : 

7 i * = S ifc -(e/x-lKu* ) (A16) 
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whose inverse is 



Hk = 9ik + (l - — )uiU k , (-417) 
allows one to write eq. (A14) in the form: 

Since g = —det(gik), we shall pose 7 = —det^ik). The ratio 7/(7 is an invariant, and its 
calculation can be performed in the co-ordinate system in which u 1 = u 2 = u 3 = 0; one 
finds: 

9 

7 = — » 
e/i 

so that eq. (A15) can be rewritten as 

H k = J^y/ll ir l ks F r8 , (-419) 



which, apart from the constant factor yj e/ fi, is just the constitutive relation for the vacuum 
case in general relativity, when 7^ acts as metric. 

We shall henceforth enclose in round brackets the indices which are either moved with 
Y k and 7^, or generated by performing a Hamiltonian derivative with respect to the latter 
tensors; therefore eq. (A19) will be rewritten as 



Due to the strict analogy with the vacuum case made evident in this way, selecting the 
form of the Lagrangian density for the electromagnetic field in the medium under question 
is for Gordon a straightforward matter. He writes: 

L' = \^^F^F lk - sVi, (-421) 



where Fik can be defined by: 

1 

7r 



F ik = <Pk,i - <Pi,k + — e (i)(k) mn (^n,m ~ ^m,n)- (-422) 



Equating to zero the independent variations of f U dS with respect to <f>i and to ipi will 
produce Maxwell's equations, a priori complemented by the constitutive relation (A14). 

It is now easy to derive the energy tensor for the electromagnetic field by starting 
from the derivation that one performs in vacuo, when the metric field 7^ is present. In 
that case Hilbert's procedure leads to write: 

5L = ±T (l)(fc) 5 7 lfc , 
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hence one gets 

T (i) (fc) = ^(F tr fWW - i^F rs FWW). {ASS) 
For the medium contemplated by Gordon one may be tempted to write: 

SL' ee iT' (l)(fc) ^, (A24) 
where L' is given by eq. (A21), and finds 

T' (i) (fc) = F ir H kr - ±5>F rs H rs , (A25) 

which is just the general relativistic version of the form proposed by Minkowski. Henceforth 
we shall drop the prime, since its omission will not lead to confusion. 

(k) 

However, can not be the energy tensor density that we are looking for, because 

it is defined with respect to the effective metric 7^, not with respect to the true metric 
gik, the only one entitled to account for the structure of space-time and, via the Einstein 
tensor, for its stress-energy-momentum tensor. In order to find the relation between T ik 
and T(j)( fc ) we simply need to express 5"f lk in terms of 5g lk . From eq. (A10) one obtains 
the variation of u l produced by the variation 6g rnn of the metric: 

Su* = \u l u m u n 5g mn = -\u l u m u n 5g mn , 
hence from the definition (A16) one gets: 

5 7 lfc = S{ g lk - (e/i - l)u l u k } = 5g lk + (e/x - l)u l u k u m u n 5g mn . (A26) 
Therefore, since TikSg lk = T^^S^y , we find immediately: 

T lk = T (i)(fc) + (e/x - l)u % u k T {m){n) u m u n . (A27) 

The mixed components of can be obtained by multiplying the left-hand side and 
the second term at the right-hand side of eq. (A27) by g kq , while the first term at the 
right-hand side is multiplied by 7 fc9 + (e/j — l)u k u q . Then 

T.« = T (i) (9) + (en - l){T (%){k) u k +u i T (mKn) u m u n }u q . (A28) 

In a co-ordinate system for which, at a given event, go- = r] ik = diag(l, 1, 1, —1) and 
u 1 = u 2 = u 3 = 0, it 4 = 1, the covariant four- vector within the curly brackets has the 
components T( Q )( 4 ), (a = 1,2,3). But under these circumstances eq. (A27) says that 
T( a )(4) = T a4 , hence in a general co-ordinate system one can write: 

T« = + (e/j - l){T lk u k + u{Y mn u m u n }u\ (A29) 

i.e., according to eq. (A25) 

T i k = F ir H kr - \5 k F rs H rs - (efi - l)Q t u k , (A30) 
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where Minkowski's "Ruh-Strahl" 

Q l = - (T> n + u l T mn u m u n ) ( A31) 

has been introduced. Since VL l Ui = 0, by substituting (A30) into (A31) one finds: 

fl 1 = F m H im - F m H m u l = u k F m (H lk u m + H km u l + H mi u k ) . ( A32) 

Equations (A30) and (A32) define the extension to general relativity of the energy tensor 
proposed by Abraham [4] , [5] for the electromagnetic field if the medium is homogeneous 
and isotropic when considered at rest. 

The generalized force density due to the electromagnetic field is eventually defined by 
the covariant divergence: 

U ee -T t k k . (A33) 
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